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The pinning of vortex lines by an array of nanoparticles embedded inside superconductors has 
become the most efficient practical way to achieve high critical currents. In this situation pinning 
occurs via trapping of the vortex-line segments and the critical current is determined by the typical 
length of the trapped segments. To verify analytical estimates and develop a quantitative descrip- 
tion of strong pinning, we numerically simulated isolated vortex lines driven through an array of 
nanoparticles. We found that the critical force grows roughly as the square root of the pin density 
and it is strongly suppressed by thermal noise. The configurations of pinned lines are strongly 
anisotropic, displacements in the drive directions are much larger than in the transverse direction. 
Moreover, we found that the roughening index for the longitudinal displacements exceeds one. This 
indicates that the local stresses in the critical region increase with the total line length and the 
elastic description breaks down in the thermodynamic limit. Thermal noise reduces the anisotropy 
of displacements in the critical regions and straightens the lines. 



I. INTRODUCTION 

The introduction of large-size nanoparticles of differ- 
ent shapes has emerged as the best practical way to im- 
prove the current performance of high-temperature su- 
perconductors. While in the first superconducting ca- 
bles the critical currents were limited by weak links, 
in the second-generation superconducting wires based 
on aligned YBa2Cu307 (YBCO) films this problem has 
been mostly resolved and critical currents are determined 
by vortex pinning. Impressive progress has been made 
to enhance critical currents in these films using both 
isotropicP^ and columnaJ^^ini inclusions. In spite of this 
progress, our understanding of strong pinning mecha- 
nisms is far from satisfactory. Theoretical estimates de- 
scribing the pinning of vortex lines by an array of strong 
pins at low temperatures were elaborated by Ovchinnikov 
and Ivlev.^^ This theory was applied to describe the be- 
havior of the critical currents in the real YBCO films 
in Refs. T2" and TS^, In particular, frequently observed 
power-law decay of the critical current as a function of 
the magnetic field with a power slightly larger than 1/2 
is naturally explained by this theory. More recently, it 
was also argued that strong pins of unknown origin de- 
termine critical currents at low magnetic fields in several 
iron pnictide compounds. It is not clear, however, to 
what extent available qualitative estimates describe the 
real situation. Due to the obvious importance of strong 
pinning by large-size inclusions for real superconducting 
materials, it is desirable to elaborate a quantitative the- 
ory describing pinning in such situations. Moreover, the 
very important issue of pinning suppression by thermal 
fluctuations does not have any theoretical description in 
the strong-pinning regime. 

In this paper we consider the pinning of vortex lines 
in a superconductor containing insulating inclusions with 
lateral sizes larger than the coherence length. We focus 




FIG. 1. Vortex line trapped by strong pinning centers at (a) 
zero and (b) finite current. 

on the pinning of individual vortex lines corresponding 
to small magnetic fields. Pinning occurs via the trapping 
of finite-size segments of a vortex lin^m with a typical 
length L, as illustrated in Fig. [l] The critical current is 
determined by the length of the trapped segment L and 
the pin-breaking force Fp, 



Therefore, in the strong-pinning regime the critical- 
current problem is mostly reduced to evaluation of the 
trapped-segment length L. In general, trapping of the 
vortex lines is a complicated dynamic process controlled 
by the competition between the pinning energy, line ten- 
sion, and intervortex interactions. Different approaches 
may be used to evaluate the trapped-segment length. 
One can assume that the pinning center always grabs a 
piece of the vortex line when it is energetically favorable. 
This assumption implies that thermal fluctuations facili- 
tate local equilibration. In this case the parameters of the 
trapped configurations can be obtained from the energy- 
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balance estimates. We will call this type of trapping the 
equilibrium regime. This energy-balance consideration 
determines trapped configurations of static vortex lines 
prepared by cooling in finite magnetic field. However, 
such consideration is not applicable to the more typical 
dynamic scenario when moving lines are trapped after the 
driving force is slowly reduced down to the critical value. 
In this case the line motion close to the critical force is 
a continuous trapping-detrapping process limited by lo- 
cal instabilities. The driven vortex can be captured when 
the line either directly collides with a strong pin or passes 
sufficiently close to it. In the second case trapping may 
occur due to the long-range pin- vortex interaction as a re- 
sult of local instability.^^ The line remains trapped until 
the force acting from the pin does not exceed the pin- 
breaking force. When the line finally stops, the pinned 
configuration is expected to be very anisotropic because 
the transverse pin-to-pin displacements which are deter- 
mined by trapping events are much smaller than the lon- 
gitudinal displacements limited by the pin-breaking cri- 
teria. 

To develop a quantitative picture and verify analyt- 
ical estimates, we explore in this paper the pinning of 
the vortex lines by nanoparticles with extensive numeri- 
cal simulations. We study the dependence of the critical 
force on the density of pins, the statistical properties of 
trapped lines including average values, and the distri- 
butions of trapping length and pin-to-pin displacements. 
We study long-range behavior of line displacements in 
the direction of the driving force and in the transverse 
direction. We also study in detail the suppression of the 
apparent critical force by thermal fiuctuations and the 
temperature dependence of the trapping parameters. 

The paper is organized as follows. In Section |TI| we de- 
scribe parameters characterizing the interaction between 
a vortex line and large-size pinning centers. In Section 
III we present analytical estimates. This includes the 
formulation of general conditions for stable trapped con- 
figurations in subsection III A[ making estimates for pa- 
rameters of trapped line in equilibrium, subsection 
and in the case of dynamic trapping, subsection 
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In the Section IIVI we describe the model used in our nu- 
merical simulations. In the Section |V| we present our 
numerical results including the zero-temperature case in 
and the finite temperature in subsection 



subsection VA 



[VBl In Section 
liminary comparisons with experiments. 



VI we discuss our results and make pre- 



with eq = ^Q/(47rAa6)^.A very important parameter is 
the pin-breaking force, the maximum force with which 
the pinning center can attract the vortex line. In con- 
trast to small defects, the pin-breaking force for large- 
size defects is limited by the line tension of the vortices. 
With increasing external force, the tips of the vortex line 
slide along the surface of the insulating inclusion until 
they meet near the equator and reconnect leading to the 
depinning of the vortex. For the in-plane current in the 
anisotropic layered material, the upper estimate for such 
a line-tension-limited force can be obtained by consider- 
ing the simple geometry of equally spaced pins aligned 
along the c axis and neglecting the interaction between 
vortex tips at the pin surface. In this case, evaluating 
the external force at which the tips meet, we obtain the 
following estimate 



Fp<(2£o/7)ln(&./en 



(3) 



where 7 is the anisotropy factor, ^min = max(^c7<5), 
is the c-axis coherence length, and s is the interlayer pe- 
riod of a layered superconductor. This force only weakly 
depends on the size and shape of the pinning center. 

The interaction of the vortex line with a remote pin is 
long ranged, due to the perturbation of the supercurrent 
ffow around the vortex by the pin. 



7r(l — ny)r 



-, for 6 < r < Aa6, 



(4) 



where Vp = {4:7r/3)bzbl is the volume of the pinning cen- 
ter and Hi are "depolarization factors" , which depend on 
the parameter jbz/b. In particular, in the case bz > b/j 
which includes close-to-spherical inclusions. 



(tanh"^ C - C) , with C ■ 
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and Ux = riy = {1 — nz)/2. 

Recently, it was demonstrated that magnetic force mi- 
croscopy can be effectively used not only for imaging but 
also for the manipulation of individual vortices^. This 
technique gives principal possibility to measure the inter- 
action between a vortex line and an individual pinning 
center and extract the relevant interaction parameters 
described in this section. 



III. ANALYTICAL ESTIMATES FOR THE 
TRAPPING OF A VORTEX LINE 



II. INTERACTION BETWEEN VORTEX LINE 
AND PINNING CENTER 

We consider first the essential parameters describing 
the interaction of vortex lines with large-size pinning cen- 
ters. The vortex pinning energy by an insulating spheroid 
inclusion with the axes b and bz is given hy^^ 

Up^2bzSoCp, Cp = \n{b/U) (2) 



A. General conditions for a static pinned line 

Consider a general vortex-line configuration trapped 
at the points (un,^n)- For simplicity, we assume that 
the forces from the pins are applied locally at the points 
z = Zn- In between the trapped points, Zn < z < 2:^+1, 
the displacement obeys the following equation 
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where / is the driving force apphed along the x axis and 
^1 = (^o/7^)>Ci is the hne tension with £i being the log- 
arithmic factor Ci = ln(rmax/^min)-'^ The displacement 
can be found as 



f{z-Zn){z-Zn^l) 

2ei 



u(2;)=Un + (Un+i-Un) 

Zn-\-l~Zn ^ti 

(6) 

The force acting from the pin on the vortex line dit z = Zn 
is given by 



F 



dw 



{Zn + 0) 



du 



{Zn - 0) 



and can be evaluated as 

p. f{Zn+l-Zn-l) 



Un+l— Un Un — Un-l 



Zn+1 ' 



-Zn-1 



(7) 



The stability condition for the trapped line is given by 



Fn < Fp for all n, 



(8) 



while the critical state corresponds to the condition that 
at least one local force reaches the pin-breaking force 



max(Fn) Fp. 



(9) 



■Top view 




One simple consequence of Eqs. ^ and ([8| is that for 
"behind" sites, Ux^n < '^a?,n-i7 '^a?,n+i7 the line-tension 
force adds to the external force meaning that they, in 
average, have shorter trapping segments z^+i — ^n, Zn — 

Zn-l- 

B. Equilibrium trapping 

Consider the trapping of a single vortex line by strong- 
pinning centers with concentration rip and pinning en- 
ergy [/p, Eq. ([2|. Assuming local equilibrium, the typical 
trapping length L and transverse displacement u are de- 
termined by the energy-balance condition^ 



FIG. 2. Upper figures illustrate a trapped vortex line in the 
metastable regime (side and top views) . The typical displace- 
ment in the direction of the force ui is much larger than 
the typical displacement in the perpendicular direction ut. 
Lower figure is a visualization of the pinned line configura- 
tion obtained in simulations (only a short section of the line 
is shown). Shadelike projections on the axis planes illustrate 
line displacements in the different directions. Short scale line 
wiggling is due to the thermal noise. The definitions of the 
trapping parameters Lt, ui^ and ut are also illustrated. 



estimate for the critical current for the equilibrium trap- 
ping 



and by the condition that the average number of impu- 
rities in the trapping volume should be of the order of 
one, 



5o . Fp 



(11) 



It is expected to increase with the pin density as -y/TLp. 



UpU L = 1. 



C. Dynamic trapping 



These equations give 



^1 ^2 
TipUp ^ 



u. 



p 



(10) 



Strictly speaking, the above conditions are obtained for 
zero current. Assuming that the trapping length does 
not change much when current is applied, we obtain an 



The equilibrium estimates for the trapping parameters 



(10) are definitely valid for the line configurations pre- 
pared by cooling at fixed field and at zero transport cur- 
rent. However, it is clear that they cannot be applied to 
the vortex lines in the critical state at low temperatures 
when moving lines are trapped after the driving force 
drops below the critical value and equilibration does not 
take place. The critical current in such a dynamic regime 
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Note that in this situation the estimate is somewhat sim- 




U 



FIG. 3. A trapped segment interacting with pinning cen- 
ter. The upper right plot illustrates the energy profile for the 
value of R below which the trapping becomes energetically fa- 
vorable but separated by the energy barrier. The lower right 
plot illustrates the energy profile near the trapping instability 
point. 



has been estimated in Ref. 11 for high fields when the in- 
tervortex interactions are essential. These estimates can 
be directly generalized to the trapping of individual vor- 
tex lines at small fields.!^ When the vortex line moves 
close to the pinning center, it may be trapped and the 
line remains trapped until the force acting from the pin 
on the vortex line does not exceed the pin-breaking force. 
In this regime two typical trapping distances, in the di- 
rection of motion, ui^ and in the transverse direction, Uf, 
are very different and have very different origins, see Fig. 
[2] These distances and the trapped-segment length L are 
connected by the geometric relation 

ripLuiUt = 1. (12) 

The longitudinal trapping distance, ui^ is determined by 
the pin-breaking condition. 



S1--F,. 



(13) 



This condition can be obtained from the x component of 
Eq. assuming that the two terms on the right-hand 
side are of the order of the pin-breaking force. 

The transverse displacement between the pins, Uf, is 
determined by the trapping events. The simplest as- 
sumption is that in most cases trapping occurs when the 
lines directly collide with the pins^^ meaning that Ut ~ b. 
This immediately gives estimates for other trapping pa- 
rameters, 



L = 



£1 



UpbFp 



(14) 



This corresponds to the following result for the critical 
current 



^0 . _ p3/2 Pp^ 

c ^"''^ ~ P V ei " 



(15) 



ilar to the result for the equilibrium case, Eq. (11), and 
has the same dependence on the pin density, oc y/n^. The 
physical assumptions behind the two estimates, however, 
are completely different. 

The above simple assumption, however, may underesti- 
mate Ut. Due to the long-range pin- vortex interaction Q 
a pinning center may capture the vortex line even without 
direct collisions. When the line passes sufficiently close 
to the pinning center it may be trapped by this center 
due to instability. To estimate the maximum transverse 
trapping distance Ut^ we consider the interaction energy 
of a segment of length L with a pinning center located at 
a distance R > b^ see Fig. H) 



E{u) = ei 



2u^ AsoVp 



{R 



with A = l/[7r(l — riy)]. This gives the interaction force 

4u 2AeoVp 

Introducing the reduced variables 



we rewrite the energy and force as 



E = ei- 



F{x) 



2R^ 
AR 



W 



(1-X)2J ' 

w 



L \^ {l-xf 

The equilibrium points are determined by 

x(l - xf = W. 

The instability point corresponds to the value of W when 
the equilibrium points vanish which happens at W > 
W^max = max[x(l — x)^] = 3^/4^. Therefore the condition 
for the instability can be written as 

AepLVp 3^ 
2eiuf ~ 44' 

which determines the maximum trapping distance Ut as 

1/4 



^[AeoLVp 
V 54£i 



(16) 



The vortex line will be trapped by the pinning centers 
located closer than this distance in the direction perpen- 
dicular to the driving force. Note that the numerical 
coefficient in this equation should not be taken too liter- 
ally because it is only correct for the simplest geometry 
illustrated in Fig. [3j Using this result, we find from Eqs. 
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(12) and (13) 



-5/4 



,1/4 



F' 



5/4 



1/4 



4/9 



4/9 



^0 > 
TipS\Fp 



1/9 



(17) 
(18) 
(19) 



This gives the fohowing estimate for the critical current 



$0 



^4/9 ^13/9 

/ton J. on 



1/9 



-5/9 



(20) 



Comparing this result with the critical current for the 
equilibrium regime (11), we can see that the two regimes 



are characterized by somewhat different dependences on 
the parameters. However, dependences on the pin den- 
sity occur to be close, the exponent in the power law, 
jc cx in this dynamic-trapping regime is somewhat 
smaller than the power 1/2 for the equilibrium regime, 
a = 4/9 0.444. 



where ap = 4^2^^ (3 ^p^/^^) and is the second 
derivative of the interaction force with respect to the line 
displacement at the maximum- force point. The veloc- 
ity criterion for the effective critical force can be written 
as Tfv = Cffc where rj is the viscosity coefficient and 
C/ <C 1. At low temperatures this gives us the following 
relation for F„ 



To 

V 



exp 



apUp 



3/2" 



Cffc 



Using the estimate fc ~ FpjLi^ the geometric relation 
( 12 ) and assuming for simplicity that Ut ~ 6, we obtain 



2/3" 



In- 



OjF^P CfToFpTlpb 



(21) 



We expect that at low temperatures the effective critical 
force and trapping parameters can be roughly evaluated 
using the simple replacement Fp Fp{T). 



IV. MODEL FOR NUMERICAL SIMULATIONS 



D. A typical pin-breaking force at finite 
temperatures 

At finite temperatures the line moves at all driving 
forces but for small forces very slow motion occurs due to 
the rare thermally activated jumps (creep regime). The 
creep and flow regimes are separated by the effective criti- 
cal force which can be evaluated using a velocity criterion. 
Such effective critical force is criterion dependent and 
thus differs from the sharply defined zero-temperature 
critical force of the depinning transition, but it has the 
advantage that it can be directly compared with exper- 
imental estimates. We do not consider creep in this pa- 
per and our purpose is to evaluate the suppression of 
this critical force by the thermal noise. The main mech- 
anism of thermal suppression is reduction of the effective 
pin-breaking force. At finite temperatures the trapped 
vortex segment has a finite lifetime on the pin even if 
the pinning force F is smaller than the maximum pin- 
breaking force. To quantify this effect, we can introduce 
the temperature-dependent force Fp{T) < Fp at which 
the trapped segments are typically released from pins. 

To evaluate this force, we assume that the line motion 
in the crossover region consists of segment jumps, mean- 
ing that the average line velocity can be estimated diS v ^ 
ui/r^ where r = tq exp [U{F)/T] is the typical time dur- 
ing which the vortex segment remains pinned and U {F) 
is the typical energy barrier for detrapping. When the 
force F acting from the pin is only slightly smaller than 
the maximum pin-breaking force, the barrier behaves as 
in the single-particle case,^^ U{F) = aFUp{l- F / Fpf^'^ , 



To develop a quantitative understanding of the strong 
pinning by the array of inclusions, we numerically simu- 
lated motion of the vortex line described by the dynamic 
equation 



eig^ + Y.F{u-Rj)5{z-zj) + e^f + FT{z,t). 



(22) 

Here / is the driving force along the x direction from the 
current, (Rj, zj) are the random pin coordinates, Ft(z, t) 
is the Langevin thermal force 



and 



F(u) 



dU{M) 



is the interacting force with a strong pin. We model the 
interaction potential by the function 



C/(u) 



52 



(23) 



where the cutoff function Gcut('^) is introduced for nu- 
merical convenience, Gcut{u) = (1 — u^/Rl^t)'^ ^ < 
Rcut and Gcut(^) = for u > Rent with R^ut > b. 
An important feature that was not taken into account 
in modeling before is the long tail in the interac- 

tion potential. However, this model does not describe 
the line-tension limited pin-breaking force. In our model 
the pin-breaking force from an isolated pin is given by 
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Fp = {3V3/S)Up/b which is achieved at = b/\^. The 
model in its original form has an unrealistic feature. For 
improbable configurations when many pins are located at 
distances smaller than 6, the vortex interaction with such 
a cluster may increase without limit. This, of course, does 
not happen in real superconductors. To bring our model 
somewhat closer to reality, we renormalized the total pin- 
vortex interaction force as Fv-p(u) = F(u — Rj) as 

p. ^-P tanh(Fv_p/Fiim) 

V — p/ lim 

SO that the maximum force cannot exceed F^^. This 
modification has only a minor infiuence on the interaction 
of the vortex line with an isolated pin. 

For numerical implementation of the model, we use the 
reduced variables 

U = OU, Z = ^^Z, t = ^^r. t. 



eib-^ au + 1 

in which the equation takes the simpler form 

5 = ^+T.^iu--R,)S{S-Zj)+eJ+FT{z,t) (24) 
i 

with 



FTAS,i)FT,a'{z',i')) = 2T5aa'5{i-i')5{z - z') 
f = T/Up 

In this dimensionless form, the equation depends only 
on the reduced temperature and reduced pin density 
= {^ib^ /Up)np. The condition of the strong pin- 
ning regime is Up < Up/eib'^ corresponding to < 1. 
For the typical parameters 7 = 5, bz = b = lOnm, and 



Up — (lOOnm) 



-3 _ 



lO^cm the reduced pin density 



can be estimated as hp ^ Upb^ / {bz^'^) ~ 10 ^. In simu- 
lations we mostly used = 1. For an isolated pin this 
gives the pin-breaking force Fp — tanh(3>/3/8) ^ 0.57. 
We also used ^cut = 50 in the cutoff function Gcut('^)- 
We consider systems of size Lx x Ly x Lz^ where Lz is 
the vortex-line length, is the size in the direction of 
line motion, and Ly the size in the transverse direction. 

One of our tasks is to compute the dependence of 
the critical force on the density of the pinning centers. 
The calculation of the steady-state critical force is not as 
trivial problem as it may appear. For a finite-size sys- 
tem there is always a metastable configuration giving the 
maximum depinning force. Assuming that such a config- 
uration is always reachable from any initial condition at 
long times, this would be the critical force we seek. In 
a large system this configuration is, however, determined 
by a rare, nontypical configuration of pins. As a conse- 
quence, the maximum depinning force slowly grows with 
the increasing the system sizes and L^, meaning that 



it is not a self- averaging parameter for a fixed line length 
Lz. To avoid this, a carefully chosen (anisotropic) ther- 
modynamic limit in all directions was proposed for d + 1 
dimensions'^. The maximum depinning force of a very 
large system is also very difficult to compute. Simulations 
of Eq. (24) at zero temperature are not very suitable for 



this purpose, because when the external force is close to 
the critical value, the vortex line traps forever in the first 
metastable state it finds. It is not clear how representa- 
tive this state is and how close the corresponding critical 
force is to the maximum value. In addition, neither the 
maximum critical force nor depinning forces for a few ac- 
cidental metastable states are very interesting quantities 
from a practical point of view. The maximum critical 
force is essentially a property of an isolated vortex line. 
A more interesting quantity is the typical pinning force 
for the finite density of the vortex lines. Indeed, even at 
low densities when vortex- vortex interactions can be ne- 
glected, we expect a typical pinning force rather than an 
extreme non-self-averaging force value to determine the 
observable critical current. 

To evaluate a typical pinning force at zero temper- 
ature, instead of the fixed-force approach, we employ 
fixed-velocity simulations using approach suggested in 
Ref. [21] Namely, we replaced the fixed external force 
/ in Eq. (24) with a slowly moving parabolic potential. 



f^K [W{t) - u^{z, t)] , Wit) = Wo + Vt (25) 

Such a potential forces the vortex line to move with 
the average velocity V. Every time the line finds a 
metastable pinned state and stops, the dragging force 
starts to increase with time until it exceeds the crit- 
ical force for this state and the line resumes motion. 
This trick allows us to explore many metastable states 
and to avoid the extreme value statistics of the sample- 
dependent critical force. The typical critical force is then 
evaluated as the average force acting on the vortex line 
in the critical configurations. 



fe = {{K [W{t)-U,iz,t)]),)t 



(26) 



where (. . .)t implies averaging over the local maxima of 
the instantaneous force. The spring constant K and the 
drag velocity V have to be taken sufficiently small so that 
they do not inffuence the calculation of the critical force. 
We typically use K - 10-^-10"^ and V = 0.001-0.002. 
The spring constants satisfy K ~ in each case, as- 
suring a proper thermodynamic limit for the critical force 
and associated critical configuration'^, and the velocities 
are small enough to assure a quasistatic stick-slip motion. 
Figure |4jleft) illustrates the typical dependences of the 
force acting on the line on the displacement of its center 
of mass for different pin densities rip and line lengths A^^ . 

We also explore the velocity-force dependences at fi- 
nite temperatures using the direct fixed-force simulations 
described by Eq. (24). Even though the term critical 



force''' is widely used in the experimental community, at 
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FIG. 4. Left: Representative dependences of the force acting 
on the line versus its center of mass location obtained using 
the fixed- velocity simulations for different pin densities Up and 
line lengths Nz (the curves are marked hy Up/Nz). The verti- 
cal segments correspond to trapped states. The typical criti- 
cal forces, /c, are obtained by averaging over the local maxima 
of these curves. Right: Examples of the displacement-time 
dependences used to evaluate average velocities at fixed tem- 
perature for the parameters shown in the plot and for different 
forces. Each force is represented by three curves correspond- 
ing to different realizations of the random potential. The line 
motion becomes more and more uneven with decreasing force. 



V. NUMERICAL RESULTS 



A. Behavior at zero temperature 
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FIG. 5. The dependence of the critical force on the pin den- 
sity. The plot contains data obtained for different system 
sizes, and the legend shows Nz / Lx / Ly. The finite-size effects 
in the critical force are weak. 



finite temperatures the concept of the critical force does 
not have an exact meaning because the velocity is fi- 
nite at all forces due to the thermal creep. Nevertheless, 
one can still introduce the characteristic force describing 
crossover between the flux-flow and flux-creep regimes 
using some average-velocity criterion. At low temper- 
atures near such critical force the line motion becomes 
very uneven, see Fig. [4]^right). It spends considerable 
time in metastable traps waiting for a strong fluctuation 
which allows it to continue motion. Such line motion 
is illustrated by animation.^ As a consequence, a proper 
averaging over such events requires huge simulation times 
and/or averaging over many realizations of the random 
potential. In addition to the critical force, we explore 
statistical parameters of trapping which allow us to un- 
derstand better the pinning mechanism. We evaluated 
the average length of a trapped segment Lt and typical 
displacements along the motion direction ui and in the 
perpendicular direction Ut. We studied the distribution 
of these parameters and their evolution with temperature 
and force. We also studied the long-range wandering of 
the line in the direction of motion and in the transverse 
direction. 



For a numerical solution the reduced equation ( 24 ) has 



to be discretized both in time and in the z coordinate. We 
typically used dt = 0.05-0.1 and dz = 1 for the discretiza- 
tion steps. To study finite-size effects, the equation was 
solved for different numbers of z-axis slices, Nz = L^/dz^ 
ranging from 512 to 2048. 
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FIG. 6. Summary of the pin-density dependences of the trap- 
ping parameters for different system sizes. As in the previous 
plot, the legend shows Nz / Lx / Ly. Upper plot: The depen- 
dences of the average trapped-segment length Lt. For com- 
parison, we show the expected trapping length extracted from 
the critical force. Lower plot: The pin-density dependences 
of the average pin-to-pin displacements in the direction of the 
force (longitudinal) and in the perpendicular direction (trans- 
verse). The longitudinal displacement, {ui)^ has a consider- 
able finite-size effect with respect to the line length Nz. 
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FIG. 7. Examples of distribution functions of trapping pa- 
rameters at one pin density and for different system sizes. 



We systematically studied the behavior of the critical 
force and the properties of trapped line configurations 
within a pin-density range spanning two orders of mag- 
nitude, from 2.7 • 10~^ to 2.7 • 10~^. Figure [s] presents 
the dependences of the critical force fc on the pin den- 
sity Up for different systems sizes Nz. We found that fc 
increases with Up according to the power law fc = foUp. 
The fit gives for the power index a value slightly smaller 
than 1/2, a ~ 0.48 and the coefficient /o ~ 0.61. In 
fact, the square-root dependence fc ~ 0.71w^ also pro- 



vides a reasonable description of the data. The power is, 
however, clearly larger than the value 0.44 suggested by 
the dynamic-trapping estimates in the case when trap- 
ping occurs due to instabilities. For the used line lengths 
Nz > 512 a noticeable finite-size effect becomes visible 
only for small pin densities Up < 3 • 10~^. 

To understand the statistical properties of trapped 
configurations, we plot in Fig. |6] the pin-density depen- 
dences of the average trapping parameters defined in Fig. 
[2j the trapping length L^, and the pin-to-pin displace- 
ments along the direction of force, ui^ and in the trans- 
verse direction, Uf. In the plot Lt{np) we show the ex- 
pected value of the typical trapped segment extracted 
from the value of the critical force Ltj = Fp/ fc. We can 
see that the real trapping segments extracted from the 



configurations indeed closely follow the expected values. 
We also see that the average trapping segments are sys- 
tematically smaller than the segments which determine 
the critical force. This is a natural behavior because one 
can expect that the critical force is determined by "weak 
spots" where the trapping segments are longer than av- 
erage along the line. The difference, however, is not very 
significant. The trapping length shows a weak but un- 
expected size effect, it slightly decreases with increasing 
total line length. We found that the product UpLtUiUt^ 
which on general ground is expected to be of the order 
of unity is, in fact, slowly decreases with Up from ~ 0.18 
to 0.14. 

From the plots of the pin-to-pin displacements, we 
can see that the pinned configurations are strongly 
anisotropic, the average displacement along the direc- 
tion of the driving force ui significantly exceeds the dis- 
placement in the transverse direction Ut. This is consis- 
tent with the dynamic-trapping picture described in Sec. 
Ill C| The difference grows with decreasing pin density 
Up. In addition, the longitudinal displacement shows a 
significant size effect, it grows with an increasing system 
size Nz. Further analysis shows that this is an indication 
of growing local stresses with increasing line length and 
suggests the destruction of the vortex lines by the pin- 
ning potential in the critical state for sufficiently large 
systems. This behavior is not anticipated by the simple 
estimates. It is interesting to note, however, that this 
size effect in ui does not lead to a significant size de- 
pendence of the critical force. On the other hand, the 
average transverse displacement Ut does not show any 
size effect. It slowly grows with decreasing rip from 1.24 
at Up = 0.0027 to 2.54 at Up = 2.7 • 10~^. As Ut remains 
comparable with the defect size, the regime in which the 
transverse trapping is determined by the long 1/r^ tail 
of the pin-vortex interaction is not quite realized. This 
explains why the power index in the force— pin-density 
dependence is larger than suggested by the metastable- 
regime estimates which assume Ut ^ 1. To obtain fur- 
ther insight into the properties of trapped lines we show 
in Fig. [T] examples of the distribution functions of trap- 
ping parameters for A fixed density Up = 2.67 • 10~^ and 
different system sizes. One can see that these distribu- 
tions are characterized by long exponential tails. There is 
a noticeable probability to find segments with very large 
ui and Lf. Note that the tails have opposite size effects 
for these parameters, the probability to find a large ui 
increases with Nz while the probability to find a large Lt 
decreases with Nz. The last trend is opposite to naive 
expectations. 

To characterize the long-range displacements of the 
line, we studied behavior of the structure factors, the 
Fourier transforms of the displacement correlation func- 
tions. 



Nz 



^u^^y{z) exp{-iqz) 
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FIG. 8. Left: Transverse and longitudinal structure factors of the line for Nz = 2048 and two densities of pins. Two regions of 
power-law q dependences are clearly observed for both components. Right: Scaled structure factors for different pin densities. 



Examples of these quantities are presented in Fig. [s] (left) 
for two pin densities. Similar to local quantities, the long- 
range displacements are also strongly anisotropic. Such 
anisotropic scaling of the displacements is a general prop- 
erty of the driven lines in the critical regime independent 
of the pinning mechanism. For both components we 
clearly observe two regions of q characterized by differ- 
ent power-law dependences 5'/^^ ex For the small- 
est q^s we found the roughness exponents Q ^ 1.14 and 
(i ^ 0.45. The value of ("t < 1 foi" the transverse direc- 
tion corresponds to the line displacements increasing as 
{[uy{z) — Uy{0)]'^) oc z'^^* at large z. On the other hand, 
the value of > 1 found for the longitudinal displace- 
ments implies that the assumed elastic approximation is 
not self-consistent in the thermodynamic limit, the aver- 
age local stress {{dux/dz)'^) increases with the line length 

Lz as In this case the longitudinal line dis- 

placements grow quadratically {[uy{z) — Uy{0)]'^) = Ciz'^ 
with the coefficient increasing with the line length as 
Ci (X L^*^^^"^^'^ This provides a natural explanation for 
the strong size dependence of the longitudinal pin-to-pin 
displacement ui in the lower plot of Fig. [Oj The negative 
size effects for the trapping lengths in the upper plot of 
the same figure can also be understood. Growing local 
stress with increasing Nz forces the line to travel longer 
distances in the longitudinal direction which increases 
the probability of finding a pin separated by a smaller 
distance in the z direction. 

The found exponents are slightly different from the val- 
ues = 1 aiid = 0.5 obtained in Ref. |23]from the ap- 
proximate functional renormalization group calculations. 
However, a similar situation was found for elastic lines in 
a plane where the numerically computed index ( = 1.25 
also exceed the predicted value (" = 1 We see 
that the transverse displacements somewhat reduce the 
exponent value in the 3D case. Since the small-g' expo- 
nents are expected to be universal, i.e., independent of 



the pinning mechanism, our results suggest that the elas- 
tic description will also break down for the weak pinning 
case. Although we have considered the linear approxi- 
mation for the elastic forces, this conclusion is expected 
to hold for the full non-quadratic energy of the deformed 
vortex line.^ 

At larger q^s we observe the regime where both com- 
ponents behave as Si^t = ^t,iQ~^ giving the short-scale 
indices = Ct = 3/2. This corresponds to displace- 
ments induced by a short-range-correlated random force, 
and such behavior is actually similar to the static Larkin 
regime for weak pinning. However, the random forces in 
our case clearly have a very different origin. 

The crossover between different regimes occurs at a 
wave vector that scales approximately as This al- 

lows us to approximately collapse the structure factors at 
different pin densities into a single curve using the scal- 
ing SM = n-'^'^'Gi^tiqn-^^^). We found ai ^ 1.32 
and at ~ 0.9. These scaled dependencies are shown on 
the right-hand side of Fig. [S] Scaling works better for 
the longitudinal structure factor. This scaling is con- 
sistent with the identification of the trapping parameter 
Ltijip) as a geometric crossover length at qLt ~ 1, be- 
tween a short-distance roughness regime with exponents 
= Ct = 3/2 and a large-distance universal regime with 
exponents ^ 0-45 and ^ 1.14. This behavior is 
again very similar to the crossover between the Larkin 
regime and the random-manifold regime established for 
the case of weak collective pinning even though disorder 
is not weak in the present case. 

We can also see from Fig. [sj^left) that both structure 
factors increase with the pin density for all wave vectors. 
This means that, in contrast to elemental pin-to-pin dis- 
placements ui^t plotted in Fig. [6j the line displacements 
at fixed z grow with increasing pin density. This behav- 
ior can be easily understood. The displacements at small 
distances are determined not only by the behavior of ui^t 
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but also by the behavior of the trapping length L^. All 
these parameters decrease with increasing pin density, 
meaning that smaller displacements occur on a smaller 
length scale. The net increase of the line displacements 
at fixed z coordinate with increasing is a consequence 
of a faster Lt{np) decrease than that of ui^ti'^ip), as can 
be seen from Fig. [6j 



B. Dynamics at finite temperatures 

In this section we consider the infiuence of thermal 
noise on the dynamic response and configurations of vor- 
tex lines interacting with strong pins. Figure [9] presents 
the temperature evolution of the velocity-force depen- 
dences and trapping parameters in the critical region for 
two very different pin densities, Up = 5.33 x 10~^ and 
Up = 1.33 X 10~^. We remind that we use the pinning 
energy of a single pin as the temperature unit. At finite 
temperatures the critical force does not have an exact 
meaning, because the line velocity is finite at all forces 
due to the thermal creep. Nevertheless, we can intro- 
duce a typical force corresponding to the crossover be- 
tween the fiow and creep regimes, similar to the voltage 
criterion widely used in experiment. We use the criterion 
V = 0.05/ for this force. The first important observation 
is that, independently of criterion, the apparent critical 
force is quite strongly suppressed by the thermal noise. 
For example, as we can see in Fig.[9j for the small pin den- 
sity Up = 5.33 X 10~^ at the temperature only 5% of the 
pinning energy, the apparent critical force is already sup- 
pressed about fourfold. As one can see from the trapping- 
parameter plots, the longitudinal pin-to-pin displacement 
ui in the critical region very rapidly decreases with the 
temperature while the transverse displacement Ut slightly 
increases with the temperature. As a consequence, at 
some temperature they become of the same order. The 
longitudinal displacement typically has a nonmonotonic 
dependence on the driving force and reaches a maximum 
at some force in the critical region (the maximum-stress 
force). The trapping length decreases with decreasing 
line velocity and its value at = 0.05/ slightly increases 
with the temperature. 

Figure [lOj^a) shows the temperature dependences of 
the apparent critical forces for a wide range of pin densi- 
ties. We can see that these dependences are quite similar. 
However, plots of the relative critical forces in Fig. 10 ^b) 



clearly show that thermal suppression weakens with in- 
creasing pin density. This is consistent with the estimate 
for the temperature renormalization of the effective pin- 
breaking force, Eq. (21 ), due to the Up dependence under 
the logarithm. For illustration, we also present the real- 
temperature scales on the top axes of the plots in Fig. 
10'a,b) computed for spherical particles with 6 = 5 nm 
and typical YBCO parameters. This scale, however, is 
very sensitive to the value of b. For example, the liq- 
uid nitrogen temperature T = 77K corresponds to the 
reduced temperature T ^ 0.041 in the plot. For par- 



ticles with b = lOnm the same real temperature would 
correspond to the much smaller value T ^ 0.011. 
As we found approximately fc oc ^fn^ at T = 0, 



ac- 
cording to the estimate of Eq. (Il5| we also expect the 

relation /c ex Fp^^ . This means that, according to 
Eq. ([21]), we expect the dependence [/c(T')//c(0)]^/^ = 
1 - {^dpTf'^ with /3^= (l/aF)ln(npo/np). This is di- 
rectly verified in Fig. 10 ^c) where we observe the approxi- 
mate linear dependencies of [/c(T)//c(0)]^/^ versus T^/^. 
Moreover, as shown in the inset, the coefficient ^^p found 
from the linear fits for different pin densities indeed has 
a logarithmic dependence on n^, ^^p ~ 1.6 ln(0.115/np). 
These observations provide justification for our assump- 
tion that the reduction of the typical pin-breaking force 
is the main source of thermal suppression of the effective 
critical force. 

We consider now the influence of thermal noise on 
the long-range behavior of the line displacements. Fig- 
ures 11 ^a,c) present the evolution of the structure factor 



with increasing temperature for the pin density Up = 
1.33 X 10~^. We can see that the slope of the small-g de- 
pendences does not change indicating that the roughen- 
ing indices and are temperature independent. How- 
ever, the coeflicient is signiflcantly reduced for the longi- 
tudinal displacements and is enlarged for the transverse 
displacements. Correspondingly, as one can see from the 
line wandering plots shown in Fig. ll^b,d), the compo- 
nents of displacements in the critical region have oppo- 
site tendencies: the longitudinal displacements decrease 
and the transverse displacements increase with increasing 
temperature. As the longitudinal displacements domi- 
nate, the lines become more straight in the critical region. 
Another important observation is that the random-force 
regime Si^t oc at large q is rapidly washed out by 
thermal noise for both components and is replaced by 
the isotropic fluctuational line wandering Si^t — T /q^ . 



VI. DISCUSSION AND COMPARISON WITH 
EXPERIMENT 

Using the square-root flt of the pin-density dependence 
of critical force in Fig. [5j we can restore the Up depen- 
dence of the critical force in real units 



1,9^3/2/^ 



(28) 



This coincides with the estimate (15) which is obtained 



assuming that the transverse trapping distance Ut is on 
the order of the pin size b. This does not contradict too 
much our numerical results, because our average values of 
Ut only slightly exceed the pin size. Note that in our sim- 
ulations the pin-breaking force Fp is flxed by interaction 
with the pin while in real superconductors for large-size 
inclusions it is determined by the in-plane line energy, 
see Eq. ([3|. Substituting this estimate, our result leads 
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FIG. 9. Temperature evolution of the velocity-force dependences (bottom plots), lengths of trapped segments (top plots, right 
axes), and pin-to-pin displacements (middle plots) for two very different pin densities, Up = 5.33 x 10"^ (left) and 1.33 x 10"^ 
(right) in the regions corresponding to crossover between flow and creep. 




FIG. 10. (a) Temperature dependences of the effective critical forces for different pin densities and line lengths, (b) Relative 
suppression of effective critical forces by thermal fluctuations for selected pi n de nsities/line lengths, (c) Plots of [fc{T) / fc{0)]'^^^ 
vs T^/"^ showing approximately linear dependences in agreement with Eq. (21 ). The inset shows the pin-density dependence of 
the coefficient /3dp with the logarithmic fit. Legends in the plots (b) and (c) imply Up/Nz- The top axes in the plots (a) and (b) 
provide the approximate real temperature scales computed for spherical particles with radius 6 = 5 nm assuming Xab = (140 
nm)/v^l - (T/Tc)2 with Tc = 90 K. 



to the following estimate for the critical force 



fc = Aceoyjupb/j, 

where, assuming < 5, we estimated Ac ^ 

bA[\n{b,/s)nn{Lt/s)]-'/\ 

For preliminary comparison with experiment we use 



results of the recent paper ^ in which approximately 
spherical (Y-Gd)2 03 particles with radii ~ 4 nm were 
embedded into the YBCO films. The typical concentra- 
tion of particles was 5 x lO^^cm"^. For estimates, we 
assume the temperature-d ependent Lon don penetration 
depth as Xab = 140nm/v'l - (T/TJ^ with = 90K 
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(a) 



0.05 q"^ 



iip=1.33xlO ,N,=2048 




— 7=0.02,^0.0041 

— 1=0.05,^0.0021 




FIG. 11. Typical evolution of the structure factors (plots 
(a) and (c)) and line displacements (plots (b) and (d)) with 
increasing temperature for one pin density. Thermal fluctu- 
ations in the critical regime reduce the line displacements in 
the direction of drive and increase the line displacements in 
the transverse direction. 



and the anisotropy 7 = 5. Such parameters correspond 
to the reduced concentration fip ~ 1.3 x 10~^. From 
the above zero-temperature result, we estimate Ac ~ 3.9 
and the critical current jc ~ 1.7 x 10^ amp/cm^ which is 
close to the low-field experimental value at 5K, 2.08 x 10^ 
amp/cm^. The reduced temperature can be evaluated as 
T ^ (T/5.3 X 10^ K) /[I - {TlTc?]- At 55 K this gives 
T ^ 0.017. According to Fig. 10 ^b), at this temperature 



thermal suppression of the critical current is expected to 
be around 50%. Taking this factor and the temperature 
dependence of the parameters into account, we expect a 
critical current jc ~ 5.2 x 10^ amp/cm^, which again 
is quite close to the experimental value of 6.22 x 10^ 
amp/cm^. We can conclude that the our strong-pinning 
result for the critical current is in reasonable agreement 
with experiment. 

Our results show that the wandering of dynamically 
pinned lines is strongly anisotropic, displacements in the 
direction of motion are much larger than displacements 
in the transverse direction. In principle, experimentally 
this can be directly demonstrated using fiux visualization 
techniques, such as a scanning SQUID or Hall probes. 
With these techniques individual vortices can only be re- 
solved at small fields. Usually statically pinned vortices 
are visualized after cooling in fixed magnetic field. Nev- 
ertheless, it should also be possible to visualize shapes of 
the dynamically-pinned vortices near the boundary of the 
Bean profile, which is formed when the magnetic field is 
applied after cooling of a superconducting sample in zero 
field. We expect the vortex field profiles to be strongly 



elongated along the direction of motion. The elongation 
length along the direction of motion is given by the lon- 
gitudinal displacement ui at a distance on the order of 
the London penetration depth. This length is expected 
to decrease with increasing temperature. 

In conclusion, we developed a quantitative description 
of individual vortex lines pinned by an array of nanopar- 
ticles. We found that the critical force grows roughly as 
the square root of the pin density. This result is con- 
sistent with qualitative estimates assuming that for our 
pin-density range trapping events mostly occur as a re- 
sult of direct collisions with pinning centers. The ap- 
parent critical force is strongly suppressed by thermal 
noise. The relative suppression reduces with increas- 
ing pin density. The configurations of pinned lines are 
strongly anisotropic, displacements in the drive direc- 
tions are much larger than those in the transverse di- 
rection. The displacement anisotropy is rapidly reduced 
by thermal noise mostly due to the rapid reduction of the 
longitudinal displacements. This leads to straightening 
of the lines in the critical region. Analyzing the behavior 
of the structure factors at small wave vectors, we found 
that the roughening index for the longitudinal displace- 
ments exceeds one. This means that the local stresses in 
the critical region increase with the line length indicat- 
ing a breakdown of the elastic description in the thermo- 
dynamic limit. The behavior of the structure factor at 
large wave vectors for both directions is typical for the 
displacement induced by random force. At finite tem- 
peratures this random-force regime is rapidly replaced 
by the thermal displacements. 
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